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PREFACE 


This  study  extends  the  work  reported  in  RAND  Memorandums  RM-2280, 

A  Correlation  of  the  Critical  Conditions  for  Homogeneous  Bare  Reactors, 
and  RM-2940-PR,  A  Method  of  Correlating  Critical  Conditions  of  Homo¬ 
geneous  Bare  Reactors  Containing  a  Resonance  Absorber,  and  shows  how 
the  effect  of  reactor  temperature  on  criticality  can  be  computed  by 
simple  hand  calculation  methods  while  retaining  the  accuracy  of  multi¬ 
group  electronic-computer  calculations. 

The  methods  presented  in  this  series  of  studies  are  intended  to 
assist  engineers  in  making  feasibility  studies  of  reactor  systems  by 
providing  a  ready  insight  into  the  effect  of  important  parameters  on 
criticality  conditions. 


SUMMARY  AND  CONCLUSIONS 


TAs  part  of  a  series  of  studies  designed  to  simplify  reactor  anal¬ 
ysis  for  use  in  engineering  evaluations,  an  analysis  is  made  of  the 
effect  of  reactor  temperature  on  criticality  conditions,  and  simple 
relations  are  obtained  for  calculating  these  effects. 

The  buckling  (and  hence  reactor  size  for  criticality)  may  be  ex¬ 
pressed  as  a  function  of  the  fission,  capture,  and  scatter  cross  sec¬ 
tions  of  the  reactor  materials  averaged  over  the  flux  spectrum.  If 
these  parameter  averages  are  known  for  one  reactor  temperature,  then 
the  simple  relations  derived  in  this  study  allow  computation  for  other 
reactor  temperatures.  A  substantial  simplification  of  the  analysis 
results  from  the  assumption  that  the  major  part  of  the  influence  of 
reactor  temperature  comes  from  the  changes  in  the  cross  sections  and 
in  the  values  of  the  flux  distribution  in  the  vicinity  of  thermal  en¬ 
ergy.  Good  agreement  is  obtained  between  the  values  computed  from  the 
equations  and  the  values  obtained  from  an  18-group  machine  calculation 
performed  for  a  series  of  reactors  covering  a  temperature  range  from 


-vii- 

ACKNOWLEDGMENT 

The  reactor-criticality  results  used  to  check  the  simplified 
method  presented  in  this  Memorandum  for  computing  the  effect  of  re¬ 
actor  temperature  on  criticality  were  provided  by  George  Safonov. 


-ix- 


CONTENTS 

PREFACE . iii 

SUMMARY  AND  CONCLUSIONS .  v 

ACKNOWLEDGMENT .  vii 

SYMBOLS.. . xi 

Section 

I .  INTRODUCTION .  I 

II.  ANALYSIS . .  2 

Criticality  Equation .  2 

Effect  of  Temperature  on  Uranium  Cross-Section  Term...  4 

Effect  of  Temperature  on  Moderator  Cross-Section  Term.  7 

Effect  of  Temperature  on  Flux  Distribution .  8 

Effect  of  Temperature  on  .  11 

Reactor  Temperature . 13 

III.  DISCUSSION  OF  RESULTS .  17 

REFERENCES . 33 


-xi- 


f. 

1 

h 

N 

P 

R 

c 

T 

Y 


2 

sme 

2  * 
smf 

<7 


SYMBOLS 


Buckling,  cm 
Neutron  energy,  ev 


“  #. 

Y  (ua  -  CT  )  yi- 
iO  fu  ail  * 

i=l  1 


Ratio  of  mean  effective  value  of  a  to  its  thermal  value, 

a 

m  * 

_  CT  , 

ami-  i 

L  a 

.  ,  amo  T 
1=1 


Fraction  of  the  neutrons  born  from  fission  into  the  i*"^1  group 
Height  of  cylinder,  cm 

Number  of  atoms  or  molecules  per  unit  volume 

Flux-spectrum  parameter  (see  Eq.  (36)  and  associated  discussion) 
Radius  of  cylinder,  cm 
Reactor  temperature,  °C 
See  Eq.  (9) 

Ratio  of  neutrons  scattered  from  epithermal  group  to  thermal  group 
per  scattering  collision  in  the  epithermal  group 

Ratio  of  neutrons  scattered  from  fast  group  to  epithermal  group 
per  scattering  collision  in  the  fast  group 

Transport  mean- free  path,  cm 

Number  of  neutrons  born  from  fission 

Macroscopic  cross  section 

Macroscopic  scatter  cross  section  of  moderator  for  epithermal  flux 
Macroscopic  scatter  cross  section  of  moderator  for  fast  flux 
Microscopic  cross  section 


Neutron  flux 
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4  Thermal  flux 
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T  Total 
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I.  '  INTRODUCTION 


This  Memorandum  continues  the  development  of  a  method  for  reac¬ 
tor-criticality  analysis  which  would  be  suitable  for  exploratory  en¬ 
gineering  evaluations.  The  effect  of  reactor  moderator  materials  on 
criticality  was  analyzed  in  Ref.  1,  and  the  effect  of  introducing  a 
resonance  absorber  was  analyzed  in  Ref.  2.  The  effect  of  reactor  tem¬ 
perature  on  criticality  is  investigated  here. 

The  method  developed  permits  the  use  of  simple  desk-calculation 
techniques  in  the  determination  of  the  effect  of  reactor  parameters 
and  reactor  temperature  on  criticality,  and  it  provides  insight  into 
the  magnitude  of  these  effects.  It  makes  use  of  a  flux-distribution 
parameter  which,  for  a  given  value,  yields  the  same  spectral-flux  dis¬ 
tribution  for  reactors  of  different  materials.  Simple  equations  are 
derived  from  an  approximate  analysis.  Their  accuracy  in  predicting 

criticality  is  checked  by  a  comparison  with  the  results  of  an  18-group 

7  235 

machine  calculation  of  a  series  of  Li  H  -  U  reactors  at  various 
moderator-to-uranium  mole  ratios,  and  at  reactor  temperatures  from 
25°C  to  1450°C.  The  data  used  in  establishing  these  correlations  con¬ 
sist  largely  of  results  of  18-  and  53-group  machine  analyses  made  by 

(3) 


George  Safonov. 
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II.  ANALYSIS 


CRITICALITY  EQUATION 

With  the  notation  of  Refs.  1  and  2,  the  equations  for  an  m-group 
solution  of  a  reactor  may  in  general  be  represented  as  m  simultaneous 
equations  of  the  form 
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4.  +  f. 

1  L 


z 

j-1 


fj  j 


i-1 

+  Z  sj. 


j-1 


si 


4 .  -  E*.  $.  -2.4.  =0 
j  si  i  ai  i 


(I) 


The  m*"*1  group  is  the  thermal  group  and  discrete  parameters  in  the  ther¬ 
mal  group  will  be  given  the  subscript  o. 

X  =  the  transport  mean  free  path  of  group  i 
4^  =  the  neutron  flux  of  group  i 

f^  =  the  fraction  of  the  neutrons  that  are  born  in  group  i 


(by  definition,  X  f.  “  1) 
i“l  1 


v  =  the  number  of  neutrons  released  per  fission  capture  of 
a  neutron 

235 

*  the  macroscopic  fission  cross  section  of  U  for  neutrons 
J  in  group  i 

S"' ,  =  the  macroscopic  cross  section  for  the  slowing  down  of 
neutrons  from  group  j  to  group  i 

•jf 

=  the  macroscopic  cross  section  for  the  slowing  down  of 
neutrons  from  group  i  to  a  lower  group;  this  quantity  has 
a  value  of  zero  for  the  thermal-energy  group 

£  .  =  the  macroscopic  capture  cross  section  for  neutrons  in 
group  i 


When  these  m  equations  (Eq .  (1))  are  added,  the  summations  of 


the  third  terms  and  fourth  terms  cancel  each  other  because  the  total 
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number  of  neutrons  slowed  into  all  energy  groups  is  equal  to  the 
total  number  of  neutrons  slowed  out  of  all  energy  groups. 

This  summation,  with  the  further  assumption  that  the  neutron - 
energy  spectrum  is  independent  of  position  in  the  reactor,  leads  to 
the  following  equation 
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$ 

r 

$ 
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I 

■i  =1 


(2a) 


where 


z 

i=l 


(2b) 


and 


X 

t 


X 

ti  $T 


(2c) 


Equation  (2a)  may  be  written  in  the  usual  form  involving  the 
buckling,  B 


V2  $  +  B2  §_  =  0 

T  T 


(3) 


where 
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i=l 
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ai 
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(4) 
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For  the  case  under  discussion,  the  reactor  is  made  up  of  U' 
a  moderator,  and  a  l/VlT  absorber ,  and  Eq.  (4)  may  be  written 
,2 
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\  ii‘  1,1 

tB  .  v 

3N  ~ 


(ua,  -  a 
fui  aui 


x  .  IN  ^ — v  *  .  IN  r v 

v  i _ m  \  i  x  \ 

i  N  ami  $  N 


i=l 


i=l 


i=l 


(5) 


where 


o  =  the  microscopic  fission  cross  section  of  U' 
=  the  microscopic  absorption  cross  section 

-  TT235 

u  =  U 

m  =  moderator 
x  =  l/v'E  absorber 


235 


EFFECT  OF  TEMPERATURE  ON  URANIUM  CROSS-SECTION  TERM 

If  the  temperature  of  the  reactor  is  changed,  the  capture  cross 

sections  in  the  thermal  group  and  the  spectral-flux  distribution  are 

changed.  The  effect  of  reactor  temperature  is  important  for  highly 

thermal  reactors  and  tends  to  decrease  as  N  /N  decreases.  At  the 

m  u 

low  values  of  N  /N  where  the  relative  number  of  captures  of  neutrons 
mu 

at  thermal  energy  is  small,  the  effect  of  reactor  temperature  on  the 
value  of  the  terms  in  Eq.  (5)  is  negligible.  Thus,  in  obtaining  ex¬ 
pressions  for  the  effect  of  reactor  temperature,  simplifying  assump¬ 
tions  that  are  accurate  in  the  highly  thermal  end  of  the  N  /N  range 

m  u 

will  be  made. 

The  first  term  on  the  right-hand  side  of  Eq.  (5)  can  be  written 
in  terms  of  its  value  at  a  reference  temperature  and  of  corrections  for 
the  effect  of  the  change  in  temperature.  A  change  in  the  percentage 
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of  thermal  flux  is  accompanied  by  an  equal  but  opposite  change  in  the 

remainder  of  the  flux.  The  change  in  the  remainder  of  the  flux  is 

spread  out  over  a  large  range  of  energy  levels  so  that  any  one  energy 

band  receives  only  a  small  contribution.  Although  there  are  resonance- 
235 

energy  bands  in  U  that  have  a  substantial  fission  and  capture  cross 
section,  these  receive  so  small  a  contribution  of  the  flux  increment 
that  they  have  only  a  small  effect  on  the  summation  term  under  dis¬ 
cussion.  Thus,  a  good  approximation  can  be  obtained  for  the  effect 
of  temperature  by  considering  only  the  contribution  of  fissions  and 
captures  in  the  vicinity  of  the  thermal-energy  level. 

Consider  a  change  in  reactor  temperature  from  a  reference  tem¬ 
perature,  T  ,  to  a  higher  temperature,.  T.  In  the  case  where  the  re- 

3L 

th. 

actor  temperature  is  T  ,  the  m  group,  which  is  the  thermal  group, 

cL 

fch 

will  be  taken  as  corresponding  to  the  temperature  T  and  the  m  -  1 

sl 

group  will  be  taken  as  having  an  energy  corresponding  to  the  higher 

temperature,  T.  For  the  reactor  at  temperature  T,  the  thermal  group 

will  correspond  to  the  same  energy  as  that  of  the  m  -  l*"^  group  for 

the  reactor  at  temperature  T  but  will  still  be  designated  the  m*"*1 

th 

group  to  retain  the  notation  that  the  m  group  is  the  thermal  group. 

If,  in  accordance  with  the  previous  discussion,  for  the  first 
term  on  the  right-hand  side  of  Eq.  (5)  we  consider  only  changes  coming 
from  groups  whose  energies  correspond  to  the  temperatures  T^  and  T,  we 


obtain 
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where  the  subscript  a  designates  parameters  associated  with  the  ref¬ 
erence  reactor  temperature.  Now  we  assume  that 


$  T 
m  -  1 


T  /a 


a  A, 


$  A  /  § 

-f* )  *AB  1 

T  /a  V  T  1  a 


(7) 


where  A  and  B  are  constants.  The  nearly  linear  relation  between 

§  ./§_  and  ®  /4>  will  be  illustrated  later.  Furthermore 

epi  T  o  T 


U°fu(m-1)  °au(m-l) 


m-1 


UCT  -  o 

fuo  auo 


(8) 


Substituting  from  Eqs .  (7)  and  (8)  into  Eq.  (6),  we  obtain 


-  a 

$ 

o 

a 

fuo 

auo 

T 

fuo  - 

auo 

r  a  -  y)  +  y 

T/a 


(9) 


f3<p)  «  Y  ■  w  r 


i=l 


where 
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and 


the  quantity  Y  may  be  conveniently  plotted  as  a  function  of  CTfuc/ (Gfu0)  • 

In  the  determination  of  Y  the  energy  interval  between  temperature  T^  and 

T  can,  if  it  is  large,  be  broken  into  more  than  two  groups.  A  plot  of 

Y  against  CTfuC/(G£u0)a  obtained  from  the  18-group  analysis  of  series  of  a 
235 

LiH  -  U  reactors  will  be  shown  later,  and  the  accuracy  of  Eq.  (9)  in 
computing  the  effect  of  reactor-temperature  change  on  the  value  of 


m 


X  <UO£ui 

i=l 


-  a 


aui 


i>  IT 

T 


will  be  illustrated.. 

EFFECT  OF  TEMPERATURE  ON  MODERATOR  CROSS-SECTION  TERM 

The  effect  of  temperature  on  the  terms  in  Eq.  (5)  for  the  1//E 
absorbers  will  be  illustrated  for  the  case  of  the  term 


X 


*. 

1 

a  ■  nr 

ami  «T 


For  a  1//E  absorber  this  term  may  be  written  as 


X 
1  =1 


i. 

i 


(10) 


If  again  we  consider  only  groups  in  the  thermal-energy  range  and  follow 
the  same  procedure  applied  in  the  case  of  the  previous  term  in  Eq.  (5), 
we  may  write,  analogous  to  Eq.  (9), 
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Using  the  definition  of  F^  given  in  Refs.  1  and  2,  namely 


E 


F4  " 


= S  K  ^ 


(ID 


we  may  write  Eq .  (11)  as 


o  F 
amo 


4 


(a  F.) 

amo  4  a 


amo 


(a  ) 
amo  a 


Y)  +  V 


(12) 


The  value  of  y  obtained  for  the  previous  case  as  a  function  of 
Cfuo/^fuo)a  aPPlies  also  in  the  case  of  E9-  (12)  if  CTfu0/(CTfuo^a  iS 


replaced  by  a  / (CT  )  , 
amo  amo  a 


An  equation  analogous  to  Eq .  (12)  can  also  be  written  for  the 
last  term  of  Eq.  (5). 

In  order  to  apply  Eqs .  (9)  and  (12)  it  is  necessary  to  determine 

w 


EFFECT  OF  TEMPERATURE  ON  FLUX  DISTRIBUTION 

An  expression  for  the  value  of  $  / relative  to  ($  / $_)  which 

o  T  o  T  a 

is  sufficiently  accurate  for  the  present  purpose  can  be  obtained  by  a 
3-group  analysis.  Again  let  it  be  assumed  that  the  major  difference 
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occurs  in  highly  thermal  reactors.  The  leakage  and  capture  of  neutrons 
in  the  fast  and  epithermal  ranges  will  be  neglected.  The  neutron  bal¬ 
ance  equations  for  the  three  energy  groups  may  be  written  as 


u  E,  4 
fuo  o 


^f  ssmf  $f  =  ° 


Mv  E  ,  4, 

f  smf  f 


y  E  4  =  0 

e  sme  e 


(13) 

(14) 


-B2X  4 

- 1—2-  +  y,  E  4  -  E  4  -  E  4=0  (15) 

3  e  sme  e  auo  o  amo  o 


where  4f  is  the  fast  flux,  4^  is  the  epithermal  flux,  and  (J^  and 
are  the  removal  fractions  for  the  fast  and  epithermal  group,  respec¬ 
tively.  The  fast  flux  is  assumed  to  cover  energies  from  10"*  to  102  ev, 
which  represents  the  range  of  energies  for  most  neutrons  at  birth  from 
fission.  The  epithermal  flux  includes  the  remainder  of  the  energy  spec¬ 
trum  between  the  fast  and  thermal  bands.  From  Eq .  (14),  we  find 


4 

e 


smf 

E 

sme 


Equation  (13)  yields 


TJ  E 


fuo 


^f  Ssmf 


From  Eqs .  (16)  and  (17),  we  have 


v  E 


fuo 


y.  E 

e  sme 


(16) 


(17) 


(18) 
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By  definition 


{  +  $  +  i. 
o  e  f 


(19) 


or 


i 

i 


X 


o 


*  *  r 

l+t  +  T 


(20) 


Substituting  from  Eqs .  (17)  and  (18),  we  find 


r-i  +  tE. 


fuo 


1  +  1 


LL  1C  u  E 

sme  smf 


(21) 


For  the  reactor  at  the  reference  temperature 


rj  • 1  +  » 

o/a 


1  .  1 
p*  2  pT*S  _ 

e  sme  f  smf 


(21a) 


From  Eqs.  (21)  and  (21a),  we  obtain 


r  - 1 

o 


fuo 


o  /a 


-  1 


(CTfuo>a 


(22) 


A  plotting  of  the  left-hand  side  of  Eq.  (22),  computed  from  machine 
multigroup  analyses,  against  (afu0)a  should  give  a  straight  line 

which  has  a  slope  of  unity. 

Thus,  if  the  results  of  a  reactor  analysis  at  a  reference  temper¬ 
ature  are  known,  then  the  criticality  conditions  at  another  tem¬ 
perature  can  be  computed  by  making  use  of  Eqs.  (22),  (12),  (9),  and  (5). 
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EFFECT  OF  TEMPERATURE  ON  \ 

The  value  of  X  in  the  thermal-reactor  range  changes  only  slowly 
with  changes  in  flux  distribution.  Thus,  a  change  in  reactor  tem¬ 
perature  should  have  only  a  minor  effect  on  X  . 

For  most  moderator  materials  the  effect  of  reactor  temperature 
on  mean  free  path  is  negligible.  For  hydrogenous  moderators  there 
is  a  sufficient  change  in  the  transport  cross  section  of  hydrogen  in 
the  thermal-energy  range  to  cause  a  small  change  in  the  transport 
mean  free  path.  This  correction,  obtained  by  considering  only  first- 
order  effects,  is  given  by  the  following  approximate  relation: 


(X  )  =  X  jr  -  (X  ) 

'  t  a  to  v  to  a 


+  X 


t  ave 


/*  \ 

$ 

1  0  ] 

o 

l  $ 

$ 

l\  T/a 

T 

(23) 


where  Xfc  i s  the  average  transport  mean  free  path  for  energy  groups 
above  thermal. 

If  we  assume  that  the  transport  mean  free  path  for  highly  thermal 
reactors  is  inversely  proportional  to  the  transport  cross  section  for 
the  moderator  c  ,  then  Eq.  (23)  may  be  written 


i  -  (X  ) 
t  v  t' 


a  <Xto>a 


(ct  )  (a  )  \  i  /  (a  )  \/$  \ 

tmo  a  _  tmo  a  j  o  j  ^  tmo  a  j/  _o  | 

CT  0  I  ”  \  1  O’  l\  / 

i  tmo  tm  ave/  T  \  tm  ave/  \  T/a 


(24) 


where 


(*  > 


to  a  N  (a  ) 
m  '  tmo  a 


1  + 


N  a 
u  tuo 

N  (CT  ) 
m  tmo  a 


-1 
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and  where  is  the  average  value  of  a  transport  crosss  section 

for  energies  above  thermal. 

7  235 

For  the  Li  H  -  U  reactors  under  consideration,  the  following 
values  for  the  transport  cross  section  of  the  moderator  were  used  in 
Eq.  (24)  to  compute  the  values  of  X  at  various  reactor  temperatures 
from  the  values  at  a  reference  temperature  of  25°C. 


Transport  Cross  Section 
of  Moderator  (barns) 

(<J  )  =  28.43 

tmo  a 


O  =9.91 

tmo 

O  =9.02 

tmo 


cr  =4.90 

tm  ave 


Temperature  (°C) 
25 
400 
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As  an  approximation  in  the  present  example,  the  following  value  was 
used : 


(X  ) 

'  to'a 


16.2  „  „ 

7a - T  =  •  Cm 

v  tmo  a 


where  16.2  corresponds  to  a  Li^H  density  of  0.82  g/cc.  This  is  based 

235 

on  the  assumption  that  the  contribution  of  the  U  is  small. 

A  comparison  of  the  values  of  X  computed  from  Eq.  (24)  with  the 
values  taken  from  the  18-group  analysis  will  be  made  later. 

Equation  (24)  does  not  include  the  effect  on  X  of  the  variation 
of  density  which  is  caused  by  change  in  temperature.  This  effect  can 
be  introduced  by  assuming  that  X  varies  inversely  with  density. 

From  Eqs .  (16)  and  (21)  we  obtain  the  following  linear  relations 
between  the  flux  components : 
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REACTOR  TEMPERATURE 

In  order  to  provide  a  consistent  method  of  relating  the  properties 
of  reactors  in  the  thermal  group  with  reactor  temperature,  the  following 
procedure  was  somewhat  arbitrarily  adopted.  It  was  assumed  that  after 
a  neutron  was  slowed  into  the  thermal-energy  group  it  would  thereafter 
be  considered  a  member  of  this  group,  regardless  of  its  subsequent  energy 
history.  Collisions  with  reactor  materials  would  then  impose  on  the 
group  of  so-called  thermal  neutrons  the  Maxwellian  energy  distribution 
which  would  include,  of  course,  the  high-energy  range.  The  values  of 
the  reactor  parameters  would  therefore  be  evaluated  on  the  basis  of  this 
Maxwellian  distribution. 

3 

The  number  of  collisions  of  type  j  per  cm  per  sec  experienced  by 
the  neutrons  in  the  thermal  group  on  the  basis  of  the  foregoing  dis¬ 
cussion  is 


$  d  E 


(28) 
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where  5  is  the  thermal  flux  per  unit  energy  at  energy  E,  and  2^  is 
the  macroscopic  cross  section  for  the  collisions  of  type  j  at  energy 


E.  The  mean  value  of  the  parameter  2.  in  the  thermal  energy  group  is 


2. 

3 


00 

I 

Q 
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EJ  dE 
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$  d  E 
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The  thermal  flux  per  unit  of  energy  on  the  assumption  of  a  Maxwellian 
distribution  is  given  by 


E  -E/KT 

2  e 

(Kir 


(30) 


For  the  case  in  which  2.  represents  a  1//1T" absorber,  let  £  be  the 

J  a 

macroscopic  cross  section  at  an  energy  equal  to  KT,  where  T  is  the 
temperature  of  the  reactor. 

Then  in  Eq.  (29) 


2 


j 


(31) 


Thus,  Eq.  (29)  may  be  written  for  the  case  in  which  2.  ■  2  is  the 

J  ® 

absorption  cross  section 
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Integration  leads  to  the  result 


E 


(33) 


Thus,  in  the  present  Memorandum,  the  mean  absorption  cross  section  for 
the  thermal  group  is  taken  to  be  equal  to  the  absorption  cross  section 
at  energy  KT,  multiplied  by  /tt/2  . 

The  thermal-energy  group  can  be  looked  upon  as  a  well  into  which 
the  neutrons  drop  at  the  end  of  their  slowing-down  process.  Once  in 
this  well,  the  neutrons  are  assigned  energies  in  accordance  with  the 
Maxwellian  distributions.  In  this  respect  the  thermal -energy  group 
is  treated  differently  from  the  higher-energy  groups.  In  higher- 
energy  groups,  the  neutrons  are  assigned  energies  within  the  energy 
bounds  of  the  group,  and  the  average  group  parameter  is  defined  by 


rii  *i  d  E 


$  d  E 


(34) 


where  2  and  4^  are  functions  of  E  and  the  integration  is  taken  over 
the  group  limit  E^  ^  to  E^. 
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The  values  of  the  thermal-fission  cross  sections  for  U  obtained 


by  means  of  Eq.  (33)  are 
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25 
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400 
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0.626 

1450 

200 

0.363 

where  ( a  )  is  the  value  of  ac  at  T  =  25°C. 
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III.  DISCUSSION  OF  RESULTS 

The  results  of  the  foregoing  analysis  will  be  compared  with  the 

results  obtained  from  an  18-group  machine  calculation  relating  to  the 

effect  of  temperature  on  the  critical  conditions  and  associated  para- 
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meters  of  a  series  of  LiH  -  U  reactors. 

The  effect  of  temperature  on  the  thermal-flux  ratio  4^/4  is 
given  by  Eq.  (22).  The  straight  line  in  Fig.  1  is  a  plotting  of  Eq. 
(22).  The  points  represent  values  taken  from  the  18-group  calcula¬ 
tion.  On  the  whole,  agreement  with  the  equation  is  good.  Some  dis¬ 
persion  occurs  at  low  values  of  N^/N^  (e.g.,  100)  where  the  assump¬ 
tions  of  this  analysis  are  inaccurate.  However,  in  this  range  the 
effect  of  temperature  on  the  parameters  of  Eq.  (5)  becomes  negligible. 
Equation  (22)  is  plotted  for  more  convenient  application  in  Fig. 

2,  where  4  / 4  is  shown  as  a  function  of  its  value  at  the  reference 
o  T 

temperature  and  of  CT£UC/  (°fuo^a ' 

Figure  3  shows  plottings  of  4^/ 4 ^  and  4^/4^,  against  4^/4^  taken 
from  the  18-group  calculations;  in  the  region  where  ^Q/^T  is  large, 
it  displays  the  linear  relation  predicted  by  Eqs .  (25)  and  (26).  Also 
as  predicted  by  Eqs.  (25)  and  (26),  reactor-temperature  variations 
have  little  influence. 

The  ratio  of  the  slopes  of  the  curves  for  4^/4,^  and  4^/4^,  gives 
a  value  of  0.82  for  iji  (see  Eq.  (27)). 

A  comparison  of  the  values  of  4^/4^,,  computed  from  the  values  of 
(VVa  by  using  Eq.  (22)  (or  Fig.  2),  with  values  taken  from  the  18- 
group  analysis  is  shown  in  Fig.  4a.  The  abscissa  in  Fig.  4a  is  the 
flux-spectrum  parameter  P,  which  in  the  present  case  is  given  by 


-18- 


Fission  cross-section  ratio,  & 

Fig.  I  —  Correlation  of  the  effect  of  reactor 
temperature  on  the  thermal -flux  ratio 
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Fig.  2— Graph  for  obtaining  thermal-flux  ratio  from  its  value  at  a 
reference  temperature  and  the  change  in  temperature  (Eg. (22)) 


Variation  of  epithermal  and  fast  flux  with  thermal  flux 
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Fig.  4a — Flux  distributions  as  a  function  of  reactor  temperature 
and  flux-spectrum  parameter  (thermal  -  flux  ratio) 
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P  =  Y  N  /N 
mm  u 
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where  Ym  =  2.6  for  a  reactor  containing  LiH  and  U  .  (A  more  de¬ 
tailed  discussion  of  this  parameter  is  given  in  Refs.  1  and  2.)  The 
reference  temperature  was  taken  as  25°C.  Good  agreement  of  the  points 

with  the  lines  is  noted  even  at  the  low  values  of  N  /N  at  which  the 

m  u 

analysis  is  invalid.  The  values  of  4^/4  and  4^/4^  are  shown  in  Figs. 

4b  and  4c,  respectively;  these  curves  also  indicate  that  the  effect 

of  change  in  reactor  temperature  tends  to  disappear  at  low  values  of 

N  /N  . 
m  u 

With  the  values  of  4Q/f^  given  by  Eq.  (22)  (or  Fig.  2)  we  can 
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now,  by  means  of  Eq.  (9),  compute  the  value  of  the  U  term  in  Eq. 

(5)  from  its  value  at  the  reference  reactor  temperature.  The  curves 

in  Fig.  5  for  temperatures  400°C  and  1450°C  were  computed  in  this 
7  235 

manner  for  the  Li  H  -  U  reactors.  The  value  of  Y  was  taken  from 
Fig.  6.  The  points  represent  values  taken  from  the  18-group  solution. 
Again  the  effect  of  reactor  temperature  tends  to  disappear  at  low 
values  of  P,  where  the  analysis  for  temperature  change  becomes  in¬ 
valid.  Good  agreement  between  the  results  of  the  18-group  machine 
calculations  and  the  present  analysis  is  obtained  over  the  entire 
range  of  values  of  P. 

The  values  of  the  moderator  term  in  Eq.  (5)  for  temperatures  of 
400°C  and  1450°C  computed  by  means  of  Eq.  (12)  and  Eq.  (22)  (or  Fig.  2) 
are  shown  in  Fig.  7.  Again  good  agreement ‘be tween  these  values  and  the 
values  taken  from  the  18- group  solution  is  noted  over  the  entire  range 
of  values  of  P.  The  values  of  F^  can  be  used  for  other  1//E  materials 
in  the  reactor,  such  as  material  x  in  Eq.  (5). 
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Fig.  4b — Flux  distributions  as  a  function  of  reactor  temperature 
and  flux- spectrum  parameter  (epithermal-flux  ratio) 
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Fig.  4c — Flux  distributions  as  a  function  of  reactor  temperature 
and  flux-spectrum  parameter  (fast-flux  ratio) 


Fig.  5  —  Effect  of  reactor  temperature  on  net  neutron -generation  term 
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Effect  of  reactor  temperature  on  moderator-absorption  term 
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Also  shown  in  Fig.  5  are  values  taken  from  a  53-group  analysis 
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of  an  HO  -  U  reactor.  In  this  case  the  value  of  Y  is  5.2.  The 
2  m 

equations  derived  in  this  Memorandum  can  be  used  as  illustrated  in 
order  to  obtain  the  temperature  corrections  for  reactors  which  con¬ 
tain  other  moderator  materials  than  those  discussed  here.  Further¬ 
more,  once  the  effect  of  temperature  on  the  terms  of  Eq.  (5)  has  been 
established,  their  values  corresponding  to  the  proper  value  of  P  can 
be  read  from  Figs.  5  and  7  for  other  reactor  materials.  This  is  dis¬ 
cussed  in  Refs.  1  and  3. 

The  values  of  used  in  Eq.  (5)  are  shown  in  Fig.  8  both  as  com¬ 
puted  by  Eqs .  (22)  and  (24)  from  the  values  of  X^  at  the  reference 
temperature  and  as  taken  from  the  18-group  solution.  The  effect  of 
change  in  density  on  Xfc  is  not  included  in  the  values  given  in  Fig.  8. 

The  application  of  the  present  method  to  the  calculation  of  re¬ 
actor  criticality  can  now  be  summarized.  Equation  (5)  may  be  written 
as 


X  B  N  N 

rrr—  =  F,(P)  -  rp  a  F,  (P)  -  rp  ct  F.  (P) 
3N  3  N  amo  4  N  axo  4 

u  u  u 


(35) 


The  flux-spectrum  parameter,  P,  can  be  computed  from 

P  =  Y  N  /N 


mm  u 


(36) 
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for  reactors  containing  only  moderator  and  U  .  Reactors  having  the 
same  value  of  P  have  nearly  the  same  flux  spectrum,  regardless  of  com¬ 
position.  (This  parameter  is  discussed  in  greater  detail  in  Refs.  1 

and  2.)  The  values  of  Y  are 

m 
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Fig. 8 — Effect  of  reactor  temperature  on  transport  mean  free  path, 

(material  density  constant  at  value  corresponding  to  25°C) 
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The  computation  of  P  for  reactors  that  include  a  resonance  absorber 
is  described  in  Ref.  3.  The  value  of  F2(P)  corresponding  to  the  ap¬ 
propriate  values  of  P  and  reactor  temperature  can  now  be  read  in  Fig. 

5.  For  reactors  having  a  value  of  P  which  is  less  than  200,  the  value 
of  F^(P)  is  independent  of  temperature  and  can  be  obtained  from  the 
curves  of  Refs.  2  and  3,  where  logarithmic  scales  are  used  for  more 
accurate  presentation  of  the  low  P-range.  The  quantity  F^(P)  corres¬ 
ponding  to  the  appropriate  values  of  P  and  reactor  temperature  can  be 
read  from  Fig.  7  and  can  be  introduced  into  Eq.  (35).  The  values  of 

O  and  cr  are  the  thermal-absorption  cross  sections  at  the  reactor 
amo  axo 

temperature.  The  value  of  X  can  be  obtained  from  Fig.  8  for  Li7H  and 
from  Refs.  2  and  3  for  other  moderators.  The  effect  of  temperature  on 
X  for  these  other  moderators  can  be  computed  from  Eq.  (24).  Thus  for 
any  reactor  composition,  i.e. ,  for  specified  values  of  N^,  N^/N^,  and 
N  /N^,  the  value  of  the  buckling,  B,  can  now  be  computed  from  Eq.  (35). 
The  effect  on  B  of  the  change  in  material  density  caused  by  temperature 
is  introduced  by  using  the  appropriate  values  of  N^,  N^,  and  N^.  The 
value  of  X  in  Fig.  8  corresponds  to  the  material  densities  at  the 
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reference  temperature;  as  the  relation  \  is  inversely  proportional 
to  densitys  this  value  can  therefore  be  corrected  to  the  density  at 
the  desired  temperature.  The  critical  size  of  the  reactor  can  be 
determined  from  the  value  of  B  by  means  of  simple  relations  for  the 
standard  reactor  shapes.  For  example 
Spherical  reactor 


Critical  Radius 


Cylindrical  reactor 


/  0-7655  \2  +  1 

U*§0  („.♦$ g2 


B 

TT 


2 

2 


where  Rc  is  the  critical  radius  and  hc  is  the  critical  height. 

The  values  in  this  study  and  those  in  Refs .  1  and  2  relate  to 
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reactors  containing  U  .  The  corresponding  values  for  reactors 
based  on  other  fissile  materials  can  be  computed  by  means  of  the 
equations  in  these  studies. 
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